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Abstract 

The Bohr-Sommerfeld quantum theory specifies the rules of quantization for circular and elliptical orbits for a one- 
electron hydrogen atom-like system. This article illustrates how a formula connecting the principal quantum number 
‘72 ’ and the length of the major axis of an elliptical orbit may be arrived at starting from the quantum mechanical 
description and how in the limit when ‘72’ is large one gets the expected classical result. 
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Introduction 


According to the non-relativistic version of quantum mechanics, the energy of the electron in a 
hydrogen-like atom depends only on the principal quantum number n (apart from the atomic 
number Z of the nucleus.) Again, in non-relativistic classical mechanics, the energy of such a 
particle, which describes an elliptical path that keeps the nucleus at one focus, depends only on 
the size (2 a) of the major axis. As one expects to recover the classical result in the large n limit 
of quantum mechanics (‘the correspondence principle’), there ought to be a connection between 
n and 2a. This connection will be examined in this article. We begin by recalling the salient 
features of the two approaches: 

The classical orbits under inverse-square type force (‘Keplerian orbits’): 

The following quantities remain conserved (Rana and Joag, 1991; Synge and Griffith, 1959): 

1. Energy, which depends only on the size of the major axis: 

*~f o) 


2. Angular momentum vector; 

3. Runge-Lenz vector, which gives the direction of the major axis. 

Hydrogen-like atom in quantum mechanics (Ghatak and Lokanathan, 1984): 

The wave function (ignoring spin) of the system : 

Vn lm (r,6,<P) = NR nl (r)Y lm (0,<P) (2) 
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where A is a normalization constant, while 

-. 1/2 


7 /m (^) = (-!)'”• 


(2/ + 1)(/ - m)! 

4 jt(l + m)! 


1 d 


l+m 


(m s> 0) (3) 


x standing for cos# ; and 


7 ;m (^^) = (-l)' ! ’7V4^^) (m< 0) (4) 


The radial part of the wave function reads 


R nl {r) = e\p l .L 2 ^(p) (5) 


L 2 J+} (p) being an associated Laguerre polynomial (Ghatak and Lokanathan, 1984) with 


p = 2 rjr 


n = 



( 6 ) 

( 7 ) 


E being the energy (E< 0) and p the reduced mass of the system. The quantum numbers n, 1, m 
have their usual meanings. 

Eigenvalue of energy ( E ) of bound states is given by 

ry2 4 

pZ e 


E = ~- 


2 2 


2 h z n 


( 8 ) 


while 


n-l-l 


y (-i) 


k +1 


[("+/)!? 


P 


{n-l-\-k)\{2l + \ + k)\k\ 


( 9 ) 


We are looking for the classical limit of the quantum mechanical solution. This requires 
that we take a large value for n. (Also large Z, which, however, is a constant for a given system.) 
But something else is also to be taken into account. In classical mechanics all three components 
of angular momentum have well-defined values. There is also a vector (Runge-Lenz vector) 
which points in the direction of the major axis; we shall come to this later. As the operators that 
correspond to the components of angular momentum (orbital angular momentum) do not 
commute with each other, only one component (z-component, by convention) of angular 
momentum and the magnitude of the angular momentum have, simultaneously, definite values in 
quantum mechanics. 

Now, writing the orbital angular momentum as Ifi, we have 

(//)+(//}+(//}-{(G)-/(/+1) (io> 
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where (^j stands for the expectation value of the operator A . Now 

= m ( - / < m < l) (11) 

with 

(/',)-(/,}-0 ( 12 ) 

If (/, \ = ±/,then (/ ; 2 ) = / 2 , and 

(?, 2 ) + (/, ! )-/ ( 13 ) 

i.e. 

(A/, ) 2 + (a/,) 2 - / (14), 

A l x , A l v being the standard deviations of l x , l v ; so that the relative uncertainty in the direction of 
the angular momentum is small. Conversely, i i'(j_ ''j = 0, then (l x 2 j + 2 |/,, 2 = /(/ + I); the 
direction of the angular momentum is now totally uncertain. Clearly the classical situation 
corresponds to (I_ ^ = ± /. For concreteness we shall take (L ''j = 1 . As a matter of fact, we shall 

consider only large values of /. Notice now that m = (I_ 'j = 1 implies 

\Y lm (0,<pf ^ (l-x 2 ) 1 ( 15 ) 

and, for large values of /, this quantity is practically equal to zero except when x = 0 or d = n 12, 
i.e.on the XY plane (Pauling and Wilson, 1935; Roy, 2010). The electron therefore spends most 
of its time on the XY plane: the first hint of the emergence of a circular or an elliptical orbit in 
the classical limit. 

Connection between ‘n’ and the length of the major axis 

We shall next find the most probable value of r for the electron; this requires finding the value of 
r (or the values of r) at which Rf tl r 2 has a local maximum (maxima.) The cases / = n — 1, n - 2, 

n - 3, n - 4 can be handled analytically. (Abel’s theorem rules out the possibility of solving 
equations of degree higher than four by algebraic techniques (Hall and Knight, 1969). 

The case l = n-1 

Here 

Rl r 1 e~ p p 2n (16) 
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and maximization of R 2 t r 2 through demanding that — (f?^r 2 ) = 0 gives p = 2n. Using this 

dp 

value of p (or, rather r = —) in the classical formula for the energy E (for a circular orbit) viz. 

2?; 

p Ze 2 
E =-, we get 

2 r 


E = - 


pZ 2 e A 

lh 2 n 2 


(17) 


this is precisely Bohr’s formula. 

Let us work out the width of the distribution of e~ p p 2 " around the maximum. Putting p : 
2n + a, and assuming that a « 2n , we have 


e p p 2n 


-2r 


(2 ri ) 


2 «g An 


(18) 


Thus the standard deviation, A p, of p is only yf2n . The relative width, ts.pl p, of the distribution 
becomes smaller with increasing n: The electron is practically confined to a circle on the XY 
plane for large n ( when m =1 and / = n — 1.) It is interesting to calculate the energy associated 
with radial motion in this case. (One does not expect any from a strictly classical stand point.) 
There are two ways of calculation: 

J{l + \)fi 2 


(i) One may subtract the energy 


2 pr~ 


linked to angular motion, from | E \ which is just the 


total kinetic energy in both classical and quantum mechanics. Setting / 
appropriate value of r viz. r = — = — we obtain 

2p ?7 


1 and using the 


(£U„=- (19) 

n 


(ii) One may, alternatively, use the uncertainty relation to first compute A p r from A r = ; and 

2?7 

then calculate ( E) nu n a i via {E) ra diai = ^ . Thus 


2 p 


2_2 


(E) 


4/T?7 


radial 


2p{Epy n 


E\ 

^ ( 20 ) 


• l=n-2 
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The primary aim of this work is to establish the connection between the size of the major axis 
(2a) and the principal quantum number (n): For a given nucleus (having charge Ze) the energy of 
the system depends only on V in classical mechanics while it depends only on the principal 
quantum number V in quantum mechanics. Remember that the major axis joins the point in the 
orbit where the radial distance of the electron from the force-centre (nucleus) is least to the point 
where the radial distance from the nucleus is most; it is also called the apsidal line. As r does not 
change for small changes in the direction of the radius vector near the major axis, the electron 
spends relatively longer periods of time at the ends of the apsidal line. Since there is no analogue 
of the Runge-Lenz vector in quantum mechanics (unless one thinks of superposition of states or 
deals with wave-packets), distances to the two ends of the major axis may be ascertained by 
finding the distances to the two maxima of R 2 , r 2 at which it has the two highest values. 



Fig. 1 


Fig. 2 



Fig. 3 
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Fig.l. For / = n— 1, R 2 , r 2 has its maximum value at r = r n = —. The distribution has a 

V 


width of approximately 



V 


Fig.2. A three-dimensional version of Fig. 1. 

Fig.3. For l < n - 2, R 2 , r 2 has two or more maxima. The values of r at which R 2 d r 2 has its 
two highest magnitudes lie at the ends of the major axis. 

In fact, these correspond to the largest and the smallest values of r at which R 2 t r 2 have 

its maxima. Note that the number of maxima of R 2 n] r 2 is (n - /). In between two successive 
maxima of R 2 t r 2 there is always a minimum of R 2 , r 2 where the function also goes to zero. 
Notice that we are relating maximum time spent to maximum of R 2 , r 2 , presumably a sensible 
guess. 


For l = n- 2, 

C 1 (P) - L 2 2 n n -2 = (2 n -2 )\[p- (In - 2)] (21) 


Thus 

Kit r 2 ~e~ p [p" -(2n-2)p"~ l f (22) 


The maxima of R 2 d r 2 are located at 


p ± = (2n — 1) ± a/ 4/7 — 3 « 2n ± 2yfn , (23) 


while there is a minimum at p = ( 2n - 2) ~ 2n. Arguments advanced above lead to a length of the 
major axis equal to + P-1 Going back to the expression for the energy in classical 

2 // ?7 

mechanics, we find that such a length of the major axis leads to the correct value of the energy. 
Thus our expectations are fulfilled for / =n - 2. 


• l=n —3 


Here 


L 


2/+1 

!+/ 


■ [(2n -3)!] 2 [- 


1 


P ‘ 


2! (2 n-5)l (2n-4)! (2//-3)!2! 


] (24) 
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Hence the roots (and minima) of R 2 , r 2 are at 

p = (2n-3)±^(2n-3) (25) 

To find where the maxima of R 2 , r 2 he, let us put p = (2 n - 2) + a. The equation satisfied by a 
is 


a 2 - (10/7 -12 )a - (12/7 -16) = 0 (26) 

Following Cardan [ Hall and Knight, loc cit] let us set a = y + z and require: yz = 
Then 


(10/7-12) 

3 


/ - 2(3* - 4) + ^4(3* - 4) 2 - (27) 

z J - 2(3 n - 4) - ^4(3* - 4) 2 - (28). 


When /7 is large 


m 

y~e 6 ( 


10/7 

3 


4) 2 (l,tu, or) 


(29a) 


and 


z ® e 


,10/7 i., 2\ 

6 (— -4) 2 (l,m, m 2 ) 


(29b) 


where (1, rw, or) are the cube roots of unity (1). 
For large /? the required solutions are 


«i - 



JR —JR 

-4) [e ¥ . + e^] = 



(30a) 
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,10/7 ' -i 

« 2 - J(— ~4)[e 6 M + e xo ] 


—JR 

6 ,. 2 - 


_ 3(— - 4) 

V^ 3 ' (30b) 


a. 


|(T- 4) 


e b co 1 + e 6 co 


= 0 (30c) 


The sum of the largest and the smallest values of p viz. (pi +/>? ) is approximately equal to 

2n 

2(2/7 - 2) or 4/7, when n is large. The length of the major axis, by our assumption, is then — and 

T] 

we are again led to the conclusion that the energy E depends on this particular quantity. 


In this last algebraically solvable case 
C--[(2„-4)!] 2 [ ' 


• l = n - 4 


P + . P 


P J 


31(2/7-7)! 21(2/7-6)1 (2/7-5)12! (2/7 - 4) S! 


] (31) 


The minima of R 2 nl r 2 are placed atp = (2/7 - 4) + a, where a is a solution of the cubic equation 


a 3 - 3(2/7 - 4)a - 2(2/7 - 4) = 0 (32) 


The required values of a can be obtained by Cardan’s method. For large n the solutions are 


JR —JR 


a ! -7(2/7-4) [e 6 . + e 6 ] = ^3(2/7-4) (33a) 


JR —JR 


a 2 ~ -yj(2n - 4) [e 6 ,co + e 6 .nr ] = -^3(2/7 - 4) 


(33b) 




■ V(2« - 4) 


e 6 co 2 + e 6 n/ 


= 0 (33c) 
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To find the values of p where R 2 nl r 2 has its maxima, let us put p = (2/7-3) + /? . Some algebraic 
manipulations lead to the following equation for/? : 

/3 4 - 6(3 n - 5 )/3 2 - 8(5/7 - 9 )J3 + (24 n 2 -126/7 +153) = 0 (34) 

Collecting the leading terms (remember that we are interested in the large n situation) we have 

ji 4 -18///? 2 - 40///? + 24/7 2 = 0 (35) 


This equation can be solved by the technique invented by Descartes [Hall and Knight, loc cit]. 
For pedagogical reasons we give a brief outline of the method. Consider the equation 


Assume 


x 4 + qx 1 + rx + s = 0 (36) 


x 4 + qx 2 + rx + s = (x 7 +kx + /)(x 2 - kx + m) (37) 

Then by equating coefficients we get 

k 6 +2qk 4 +(q 2 -4s)k 2 -r 2 =0 (38) 

This is a cubic equation in k 2 . After solving this equation we can find the values of m and /. 
Finally x (four possible values) is obtained from the pair of equations 

x 2 +kx + l = 0 (39) 


and 

x 2 -kx + m = 0. (40) 


Applying this technique to our equation for /?, we get 


V28n ± , 28/7 - 4(5/7 - 19 ^) - V28^ ± , 28/7 - 4(5/7 + 19 ^) 

-I- s -I-^ ( 4i) 
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In the large n situation the highest and the lowest values of /:> are respectively ^-(V28 + V8)Vn 

and - ^-(-s/28 + V8)Vn . Clearly the sum of the corresponding values of p is approximately equal 
to An, as anticipated. Our stand is therefore vindicated. 






69 






















European Journal of Physics (Science) Education Vol.2 No.2 2011 ISSN 1309 7202 Syam and Roy 




European 


Fig.4. Plots (schematic) of the variation of Rf tl r 2 with p for different values of /. 

Conclusions 

The cases for other 7 ’values may be handled by numerical methods. However, the trend is 
evident, for it has to be consistent with Bohr’s correspondence principle. We hope that through 
this article we have been able to illustrate how quantum effects disappear and classical results 
emerge in the limit of large 7z’ values. It may also help the students to accept quantum 
mechanics as a deeper and refined description of nature. 
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